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ABSTRACT 

We present a process that accounts for the steep decline and plateau phase of the Swift X-Ray Telescope (XRT) 
light curves, vexing features of gamma-ray burst (GRB) phenomenology. This process is an integral part of the 
“supercritical pile” GRB model, proposed a few years ago to account for the conversion of the GRB kinetic energy 
into radiation with a spectral peak at E p y ~ m e c 2 . We compute the evolution of the relativistic blast wave (RBW) 
Lorentz factor T to show that the radiation-reaction force due to the GRB emission can produce an abrupt, small 
(~25%) decrease in T at a radius that is smaller (depending on conditions) than the deceleration radius R D . Because 
of this reduction, the kinematic criticality criterion of the “supercritical pile” is no longer fulfilled. Transfer of the 
proton energy into electrons ceases and the GRB enters abruptly the afterglow phase at a luminosity smaller by 
~m p / m,, than that of the prompt emission. If the radius at which this slow-down occurs is significantly smaller than 
Rd, the RBW internal energy continues to drive the RBW expansion at a constant (new) T and its X-ray luminosity 
remains constant until Rd is reached, at which point it resumes its more conventional decay, thereby completing 
the “unexpected” XRT light curve phase. If this transition occurs at R ~ R D , the steep decline is followed by a flux 
decrease instead of a “plateau,” consistent with the conventional afterglow declines. Besides providing an account 
of these peculiarities, the model suggests that the afterglow phase may in fact begin before the RBW reaches 
R ~ R d , thus providing novel insights into GRB phenomenology. 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are extremely bright explosions 

at cosmological distances (Costa et al. 1997; van Paradijs 

et al. 1997) with isotropic luminosities occasionally exceeding 
~10 54 erg s _1 . Their durations are in the range ~0. 1-1000 s 
and their luminosities peak at an energy close to the electron 
rest mass energy, E^ ~ 1 MeV (however, the accumulation 
of observational data has shown that this characteristic energy 
exhibits a wider distribution, ranging from as low as a few keV 
(Campana et al. 2006) to as high as 15 MeV (Axelsson et al. 
2012), in correlation with either the isotropic energy released in 
the burst, E lso , or its peak isotropic luminosity L pjs0 ). 

They are believed to originate in the collapse of stellar cores 
(long GRBs) or the mergers of neutron stars (short GRBs), 
processes that result in jet-like relativistic outflows with Lorentz 
factors T < 300 (but on occasion exceeding values 500-1000; 
Abdo et al. 2009a, 2009b; Ackermann et al. 2010; Hascoet et al. 
2012). It is generally considered that the kinetic energy of these 
outflows is converted efficiently into radiation in collisions of 
shells of matter ejected at different times by the GRB “central 
engine”; such collisions are thought to be necessary in order 
to produce the observed rapid GRB variability (see, however, 
Narayan & Kumar 2009) and spectra with characteristic GRB 
signatures, i.e., a maximum luminosity at energies E pk ~ m e c 2 
(for reviews, see Piran 2004 for bursts prior to the launch of Swift 
and Zhang 2007 for bursts after the Swift launch). Following 
this most luminous, prompt, y- ray emission phase, GRBs shift 
into their afterglow phase. In this phase, their luminosity is 

substantially lower and their peak emission shifts into the X-ray 
band. The longer duration of this phase (Ar ~ 10 5 s) allows a 


more precise localization and optical detection, which can then 
provide the GRB redshift. 

According to prevailing theory (Piran 2004; Zhang 2007), 
GRB emission is due to synchrotron radiation by electrons ac- 
celerated, in the prompt phase, in the shocks of the colliding 
shells. In the afterglow, the emission is due to electrons accel- 
erated in the forward shock of the expanding relativistic blast 
wave (RBW). As the RBW expands, it sweeps more matter. 
After it has swept up an amount M — E/c 2 T 2 (E is its total 
injected energy, T its asymptotic Lorentz factor) at a distance 
Rd, its Lorentz factor decreases and it is thought to enter the 
afterglow phase, as surmised by its declining X-ray and optical 
fluxes. Under these assumptions, one can calculate the expected 
X-ray flux decrease with time, which turns out to be a power law, 
Fx oc t~ a , with a ~ 1 in spherical (Sari et al. 1998) and a ~ 2 
in jet-like (Sari et al. 1999) flows. Indeed, the early, sparsely 
sampled, pr e-Swift light curves appeared consistent with such 
a behavior. However, more densely sampled X-ray light curves 
with the X-Ray Telescope (XRT) aboard Swift uncovered sig- 
nificant deviations from this behavior. So, following the prompt 
Swift Burst Alert Telescope (BAT) y-ray emission, typical XRT 
afterglows (Tagliaferri et al. 2005; Nousek et al. 2006; Evans 
et al. 2009) comprise a segment of much steeper flux decline 
(oc r -3 to r~ 6 ), followed by either a less steep power law (Liang 
et al. 2007), or a 10 2 — 10 5 s period of nearly constant flux (a 
“plateau”), followed finally at t = 7j, r k by the more conven- 
tional power-law decline ~f _1 . In addition. Swift follow-ups 
also discovered (Burrows et al. 2005) occasional flares on top of 
these light curves, as late as ~10 5 s after the BAT trigger. These 
unexpected details in the GRB afterglow light curves were added 
to the other already open problems related to the GRB prompt 
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emission, namely the nature of their “inner engine,” the non- 
dissipative transport of their energy to the emission region to 
distances R ~ 10 16 -10 17 cm and its efficient dissipation there, 
and the physics behind the distribution of the GRB peak energy 
£p k ~ 1 MeV (Mallozzi et al. 1995; Preece et al. 2000). 

These prompt GRB emission issues are usually settled by fiat 
in most studies, e.g., by assuming that a large fraction (~50%) 
of the proton post-shock energy is converted into electrons 
with a minimum Lorentz factor y mm , chosen so such that the 
burst luminosity maximum would appear at E p ^. However, this 
exhausts essentially most of the models’ freedom, setting the 
afterglow evolution on the path described in Sari et al. (1998), 
thus making an account of the observed afterglow light curve 
shapes all the more pressing. There have been a number of 
attempts to account for at least some of these features. Thus, 
Kazanas et al. (2007) suggested that if the post-shock proton 
distribution function comprises, in addition to a relativistic 
Maxwellian of T ~ Tm p c 2 , also a power law that extends to 
energies ym p Tm p , one could account for the steep decline 
followed by the more conventional oc f _1 power-law decline 
light curves of GRB afterglows. A similar proposal was put 
forward by Giannios & Spitkovsky (2009), who employed the 
realistic electron distributions produced in particle-in-cell (PIC) 
simulations, which do have a form similar to that conjectured in 
Kazanas et al. (2007). More recently, Petropoulou et al. (2011) 
provided models with the desired general afterglow shape by 
adjusting the maximum electron distribution Lorentz factor y max 
in such a way that the steep decline represents synchrotron 
emission by electrons near y max (which are cooling fast), 
while attributing the constant X-ray flux component to inverse 
Compton emission by the (much more slowly varying) lower 
energy section of the electron distribution. Another commonly 
accepted interpretation of the steep decay phase in the GRB early 
afterglow light curves is that the high-latitude emission, i.e., the 
curvature effect (Zhang 2006 and references therein), follows 
the prompt emission phase. The goal of the present work is to 
show that the XRT afterglow light curves can be accounted for 
by incorporating the supercritical pile model (hereafter SPM) 
and its feedback on the dynamics of the GRB RBWs that gives 
rise to the GRB. 

The SPM (Kazanas et al. 2002; Mastichiadis & Kazanas 2006, 
2009) has been introduced to address the issue of the GRB dis- 
sipation and the apparently efficient conversion of the RBW 
kinetic energy into radiation with the observed spectral charac- 
teristics. The fundamental process of this model is a radiative 
instability that can convert the internal energy of the RBW rel- 
ativistic protons into relativistic e + , e~ pairs. The conversion 
takes place on timescales A R/c (A R is the typical width of 
the RBW) via the py -> p e + e~ reaction, provided that certain 
kinematic and dynamic threshold (criticality) conditions are ful- 
filled, which are discussed in the next section. Unlike the more 
conventional GRB models in the literature, the SPM does not 
require (but does not forbid) accelerated particle populations be- 
sides those produced by the isotropization of the RBW kinetic 
energy behind the shock. Most importantly, it does not invoke 
an ad hoc equipartition between the proton and electron energy 
densities behind the shock, as it includes the dynamics that con- 
vert the proton energy into e + , e ~ . Finally, a consequence of the 
kinematic threshold condition of the p y — > p e + e~ reaction 
is the natural emergence (after all relevant Lorentz transforma- 
tions) of a characteristic photon peak energy, Zs pk ~ 1 MeV 
(in the observer frame, assuming that the process operates 
close to its kinematic threshold), largely independent of the 


RBW Lorentz factor and in agreement with observations (e.g., 
Mallozzi et al. 1995; Goldstein et al. 2012). In more conven- 
tional models, such a characteristic photon energy occurs only 
at the expense of assuming the presence of a low energy cut-off 
in the electron distribution function. 

More recently, Mastichiadis & Kazanas (2009) explored 
numerically the SPM from the prompt into the afterglow stage 
by computing the evolution and spectrum of a RBW of Lorentz 
factor F o = 100 propagating through a medium of density 
n(r) oc r~ 2 , representing the wind of a Wolf-Rayet star. In 
this treatment, they incorporated the effects of the radiative 
drag into the RBW evolution equations, where these effects 
are introduced in the production of the GRB spectra by the 
bulk Comptonization of upstream scattered photons, a crucial 
element of this model. They showed that when both threshold 
conditions of the model were fulfilled and the energy stored 
in protons is converted into radiation, the resulting radiation 
reaction force reduced the RBW Lorentz factor T and also the 
GRB flux. In the particular case they examined, the drop in T 
was sufficiently large to render the RBW subcritical. The GRB 
thus entered the afterglow phase (with its luminosity coming 
only from the electrons being swept up by the RBW) after only 
a couple of seconds (in the observer’s frame), thus producing 
a short GRB, even though the RBW was assumed to propagate 
through a pre-supernova stellar wind medium. 

In this work, we employ a simplified version of the RBW 
evolution discussed in Mastichiadis & Kazanas (2009) in order 
to focus our attention to the entire evolution of T from its 
acceleration phase into its late time decline. Using this simplified 
version, we provide an account of the vexing steep decline and 
plateau phase in the early GRB afterglow light curves observed 
in most long GRBs. The simplified version of this approach 
allows us to compute the evolution of both conical and parabolic 
GRB jet configurations in a medium of constant density p from 
their accelerating phase (i.e., T ~ 1) to their adiabatic decay past 
the deceleration radius. In Section 2, we introduce the general 
framework of the SPM and criticality conditions. In Section 3, 
we present our results and then in Section 4 we summarize our 
findings and present our conclusions. 

2. THE SUPERCRITICAL PILE MODEL IN BRIEF 

The process described in this section was first used in 
the context of active galactic nuclei (AGNs) by Kazanas & 
Mastichiadis ( 1 999), who employed it to argue for the possibility 
of a hadronic origin of the relativistic electrons in blazars. It 
involves the combination of: (1) The relativistic proton plasma 
radiative instability of Kirk & Mastichiadis (1992) as applied 
to a RBW. (2) The increase in the energy of synchrotron 
photons produced in this plasma, which, upon their scattering 
in upstream-located matter (referred to as the “mirror”), are 
then re-intercepted by the RBW, as discussed in Ghisellini & 
Madau (1996). The instability of Kirk & Mastichiadis (1992) is 
basically that of a nuclear pile: synchrotron photons, produced 
by e + , e~ pairs, interact with the relativistic protons of the 
plasma to produce more e + , e~ pairs; the set-up is radiatively 
unstable if the column of the plasma is sufficiently large such 
that at least one of the A photons ( N ~ y/by 2 = 1 /by, y is 
the electron Lorentz factor and b the magnetic field normalized 
to the critical one B c — m 2 c 3 /eh — 4.4 x 10 u G) produced by 
a member of an e + , e~ pair of energy E e ~ ym e c 2 produces 
another pair in reaction with a relativistic proton before escaping 
the system by the process py —> p e~ e + . 
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If R is the size of the plasma and n f} is the proton density, 
this last constraint reads a py Rno > 1 / N ~ by. However, the 
py process requires that the energy of the synchrotron photon 
E s be sufficiently high to produce a pair in the proton rest 
frame, i.e., y E s > 2m e c 2 . Considering that E s ~ by 2 m e c 2 , 
the kinematic threshold reads y 3 > 2/b. Incorporating this 
in the column density constraint, one obtains o Rn lt y 2 > 2. 
Applying these considerations to the particles in the post-shock 
region of a RBW, one can set their y values equal to the RBW 
Lorentz factor T so that the criticality conditions are expressed 
in terms of the bulk RBW kinematic properties, i.e., F > 2/b 
and (i pY Rri()Y 2 > 2. 

If the synchrotron photons scatter upstream of the RBW in a 
“mirror” (in AGNs, these mirrors are broad-line region clouds) 
of scattering depth r nmT , upon their re-interception by the RBW 
they have energies larger by 4r 2 . This modifies the kinematics 
and also the column density conditions to 

T 5 > (1/2*) and 2T mirr n 0 a RF* > 1. (1) 


With the above setting for the conversion of proton energy into 
electrons and photons, it was shown in Kazanas et al. (2002) 
that the threshold of the py — > p e~ e + reaction translates into 
an energy bF 5 in the observer’s frame, which, by the first of the 
relations above, implies a peak emission energy at roughly the 
electron rest mass, in agreement with observations. 

The observed peak energy occurs at this value only if the 
process operates close to the kinematic threshold at all times. 
The fact that bF 5 ~ 1/2 does not suffice to produce a burst be- 
cause, while bF' may be well above the threshold, say 10 MeV 
or higher, rapid proton energy release requires also that the 
column of swept-up protons be sufficiently high such that the 
dynamic threshold in Equation (1) also be satisfied. Apparently, 
the larger for'' is, the faster the accumulated energy will be 
released once supercritical. Therefore, the model does not ex- 
clude the higher values of Tjy observed recently (Axelsson et al. 
2012; Guiriec et al. 2013). More importantly, the bright bursts 
analyzed in these references indicate a correlation between the 
value of £pt and the burst luminosity (L oc £“ k , a ~ 1.3) during 
the evolution of the same burst (not addressed to the best of our 
knowledge by models to date), suggesting additional nuances; 
it would be unwise to speculate at this point. 

On the other hand, values of T pk smaller than 1 MeV 
(Campana et al. 2006) averaged over the burst duration may 
be due to larger viewing angles of the GRB jet 6 > 1/T, as 
discussed in Ioka & Nakamura (2001). Also, E p t < 1 can be 
obtained for 0 < 1/T if the shock produces, in addition to 
protons with energy ~r m p c 2 , a power-law tail that extends to 
energies y y\ > F such that by 3 > 1 /2, while £>T 5 < 1/2, 
as discussed in Kazanas et al. (2007). Independently of the 
specifics discussed above, the important point to bear in mind 
is the emergence of a characteristic energy in agreement with 
observations after all Lorentz transformations have taken place, 
as a result of the physics of the dissipation process. 

The evolution of T of a RBW is given by the coupled mass 
and energy-momentum conservation laws (see, e.g., Chiang 
& Dermer 1999). In the case that the RBW plows through a 
radiation field, one must also include the effects of the radiation 
reaction of the RBW as it plows through the radiation that has 
scattered upstream of the shock (Mastichiadis & Kazanas 2008; 
Boettcher & Principe 2009); these effects are given below: 


dM 

~dR 


= 47rR 2 r p — 


E 

c^F 


( 2 ) 


and 

dF _ ^ AR 0 /R 2 _ 4jrR 2 pr 2 _ F md 
dR ~ 1 + ARu/R M Me 2 ’ 

where R(t) is the radius measured from the center of the original 
explosion and A = Eq/ Mqc 2 , where Eq and Mq represent the 
initial total energy and rest mass, respectively, of the flow at 
R = /\*o, corresponding to the radius of the GRB progenitor. 
Here, E represents the radiation emission rate as measured in 
the comoving frame and F md is the radiation reaction force 
exerted on the RBW by the radiation field exterior to the flow, 
given by 

647T , • 

F rad = T b n e o T RY A E, (4) 

9c 

where t/, is the RBW Thomson depth, ay is the Thomson 
cross section, and n e is the circumstellar matter electron density 
assumed to be the same as the proton density n d used in 
Equation (1). 

Mastichiadis & Kazanas (2009) applied these equations to 
compute both the evolution of the Lorentz factor F and the 
emitted radiation for a RBW propagating through the wind of a 
Wolf-Rayet star, i.e.,amedium with density profilen(R) oc R 2 
that presumably was the progenitor of the supernova that gave 
rise to the RBW. The values of E and F md were computed 
by implementing the numerical code, originally described in 
Mastichiadis & Kazanas (2006), to solve the equations 

3 rij 

yr 1 +L i + Qi = 0, (5) 

at 

where the functions n, represent the differential number densities 
of protons, electrons, and photons with the index i taking any of 
the subscripts p, e or y, while L, denotes the losses and escape 
and Qi denotes the injection and source terms in the system. 

The detailed calculations of Mastichiadis & Kazanas (2009), 
using the parameters R a = 10 14 cm, « 0 = 8 x 10 8 cm -3 , Tq = 
100, Bo — 4.4 x 10 4 G, and Eq = 10 54 erg different from those 
used herein — except for Eq — but appropriate for the setting 
considered, confirmed the qualitative estimates concerning the 
positions and effects of the problem thresholds. In addition, they 
showed that the radiation reaction effects were significant, and 
had an immediate effect on T, which slowed down over a length 
scale short compared with Rq to a value lower than that required 
by the kinematic threshold of the problem. This resulted in the 
precipitous decline of the GRB luminosity as the only available 
energy to be radiated from that point on was that of the swept- 
up electrons. At the same time, they computed the resulting 
spectrum and found it to be in agreement with the basic premises 
of the SPM (i.e., exhibiting a peak at E p ^ ~ m e c 2 ). The resulting 
spectrum also agreed with the spectrum softening significantly 
with the decrease in T effected by the radiation reaction process, 
with the GRB thus entering the afterglow stage. 

3. RESULTS 

In the present paper, we study a simplified version of the 
evolution of the RBW Lorentz factor F in a medium of constant 
number density «o within the SPM. To this end, we begin our 
computations at the radius of the GRB progenitor, Rq, where 
we set To = 1, with approximate estimates for the resulting 
luminosity, a fact that allows a much broader search in parameter 
space and exploration of the evolution over longer time scales. 
The evolution of F is followed from its initial accelerating phase 
to its saturation (constant T) and slow-down stages, attributing 
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each to the prompt or afterglow stage depending on whether the 
criticality conditions are fulfilled. 

To reduce the number of free parameters, we assume that the 
magnetic field is in equipartition with the post-shock pressure so 
that B — 0. 4 ('no/ 1 cnr’) l/2 r G. Then, the long-term evolution 
of the Lorentz factor depends on the free parameters Eq,Rq, and 
no, which determine the radius where the kinematic and dynamic 
conditions are satisfied so that the RBW becomes supercritical. 
It could happen that for certain parameter combinations the 
threshold conditions are satisfied at more than one radius, in 
which case the released energy should be proportional to the 
time between such bursts, as found by Ramirez-Ruiz & Merloni 
(2001). We examine both conical and parabolic configurations 
of the GRB jet with parameters Rq = 10 n cm,no= 100cm -3 , 
Mqc 2 = 5 x 10 51 erg, and total isotropic energy Eq = 10 54 erg. 

3.1. Conical Outflows 

The evolution of T is given by the solution of the coupled 
Equations (2) and (3). To simplify our treatment, instead of 
using the numerical code employed in Mastichiadis & Kazanas 
(2006) and Mastichiadis & Kazanas (2009) to calculate the 
radiation emission rate E, we use the fact that once the criticality 
conditions are satisfied almost all the energy in the swept-up 
protons is immediately radiated away, so that 

E = £; n j = 47rR 2 p(r 2 — T)c 3 , (6) 

where E mi is the proton energy injection rate (Blandford & 
McKee 1976). Letting x = R/Rq, the evolution equations in 
Equations (2) and (3) become 

= AtzR1x 2 Tp — 4jzx-rIp(T — 1) (7) 

dx 

and 

dC _ TA 4jrx 2 R 3 pT 2 

dx x{x + A) M 

2567T 2 , , . 

- T b n e cr T x 3 R 4 0 pr\r 2 - T). (8) 

As long as the threshold conditions in Equation (1) are not satis- 
fied, then the evolution of the RBW is described approximately 
by Equations (7) and (8) without the E and F m d terms on the 
right-hand side. This is the standard non-radiative case men- 
tioned earlier in which the Lorentz factor reaches its asymptotic 
value r ~ A = Eq/Mqc 2 = 200 and then proceeds with the 
conventional decline of afterglow theory, shown by the dashed 
curve in Figure 1 . For the chosen values of the GRB parame- 
ters, the kinematic condition is satisfied at log(R/Ro) — 2.77 
(represented by the first vertical red line in Figure 1) before the 
Lorentz factor reaches even its asymptotic value. Once enough 
matter is piled up such that the column of accumulated hot pro- 
tons exceeds its critical value, the RBW becomes supercritical 
at \og(R/Ro) — 4.55 (represented by the green vertical line 
in Figure 1); the proton accumulated energy is released on the 
shock light crossing time scale and results in a sudden drop in 
r of the RBW due to radiative drag. The decrease in T reduces 
the value of bT 5 below its threshold value at log( R j Rq) ~ 4.56 
(represented by the second red line in Figure 1) and arrests 
the conversion of proton energy into radiation. The luminos- 
ity drops precipitously, by roughly a factor m p /m e (m p and m e 
are the proton and electron masses, respectively) as the emit- 
ted radiation now comes only from the cooling of the electrons 



Figure 1. Lorentz factor for a RBW propagating in a constant density 
environment for a conical GRB outflow. The dashed curve represents the 
evolution with radius in the adiabatic case without radiative drag. The red 
vertical lines represent the region where the kinematic condition is satisfied and 
the green line indicates the radius at which the dynamic condition starts to be 
satisfied. The GRB parameters are given in the text. 

(A color version of this figure is available in the online journal.) 

swept by the RBW and the GRB enters the afterglow stage (one 
should note here that, in contrast with most models, the SPM 
provides natural, physical grounds for the separation of GRBs 
in the prompt and afterglow stages). 

However, despite the decrease in luminosity that follows the 
reduction in T due to the effects of radiation reaction, the 
rest mass accumulated to this point may be too small, for 
the given RBW internal energy Eq, to produce a decrease in 
r in the manner expected generally for distances R > R « . 
Therefore, the RBW evolution must continue at a constant (or 
even increasing) F, even though, according to the premises of the 
SPM, the GRB has entered the afterglow stage (the kinematic 
threshold condition is not satisfied; the py — > p e + e~ reaction 
does not take place and hence E pk <<C m e c 2 ). The value of T will 
remain constant at a distance equal to the deceleration radius Ro, 
as shown in Figure 1 . We contend that this stage is responsible 
for the plateau phase observed in the XRT light curves. 

In Figure 3 (left) and (right), we plot the evolution of two 
RBWs with the same initial conditions as those of Figures 1 
and 2, except for the value of Mqc 2 , which now determines the 
asymptotic value of Too in the absence of radiation reaction. The 
value of Mqc 2 for these features is set to Mqc 2 — 3.5 x 10 51 erg 
and M 0 c 2 = 2 x 10 51 erg for Figure 3 (left) and (right), 
respectively, implying asymptotic T values of ~286 and 500. 
For these conditions, the radiation reaction feedback is much 
larger, with the sharp transition taking place at a distance 
increasingly closer to Rq with increasing value of Too. This then 
implies a concomitant decrease in the length of the constant T 
section in the afterglow stage, the latter effectively disappearing 
for the largest value of Too. This is of interest because of the 
correlation between the plateau luminosity and the time 7j,k of 
its break to the more conventional afterglow decrease with time 
(Dainotti et al. 2008, 2010). It is worth noting that the evolution 
of T in Figure 3 (right) would likely correspond to one of the 
typical shapes of the afterglow curves indicated in Willingale 
& O’Brien (2007), i.e., that of a steep decline, followed by a 
conventional time decrease of the XRT flux. 

3.2. Parabolic Outflows 

In a parabolic outflow, the accumulated number of ambient 
particles increases as N{R) oc R 2 , unlike the conical case 
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Figure 2. Lorentz factor for a RBW propagating in a constant density 
environment for a parabolic GRB outflow. The dashed curve represents the 
evolution with radius in the adiabatic case without radiative drag. The red 
vertical lines represent the region where the kinematic condition is satisfied and 
the green line indicates the radius at which the dynamic condition starts to be 
satisfied. The GRB parameters are given in the text. 

(A color version of this figure is available in the online journal.) 

where N(R) oc R\ Therefore, the evolution equations with 
the previous assumption that E = E ln] are given by 

dM , , 

= Ait RlxY p — AtzxRqP(T — 1) (9) 

dx 

and 

dr 2r A 47rx/?gpr 2 
dx 3x(x 2 / 3 + A) M 

- 15 ^T h n e a T x 2 R*pY\Y 2 - V). ( 10 ) 

Solving these differential equations numerically in a constant 
density medium with the same GRB parameters used previously 
in the conical case gives the evolution of the Lorentz factor T 
with radius R, as shown Figure 2. Again, the dotted curve repre- 
sents the evolution in the adiabatic case without radiative drag. 
In the parabolic configuration, the expanding RBW accumu- 
lates mass at a lower rate and therefore has a larger deceleration 


radius R D than the conical case; for the same reason, the lumi- 
nosity released is smaller. In fact, both threshold conditions are 
satisfied during the acceleration phase of the RBW when the 
first term on the right-hand side of Equation (10) is dominant, 
i.e., before the RBW Lorentz factor reaches its asymptotic value 
r = 200. These factors explain the delay and smoothness of the 
drop in T when the RBW becomes supercritical. 

As can also be seen in Figure 2, the slowing down of the 
RBW due to radiative drag during the supercritical phase reduces 
for 3 below its threshold value, thereby ending the prompt GRB 
phase. This is followed by a period of constant T, which is 
longer than that obtained previously due to the fact that the 
RBW becomes supercritical during its accelerated expansion 
phase and the fact that, for a given Eq, it sweeps up matter at 
a lower rate than a conical RBW. After all the non-adiabatic 
effects have died down and once the RBW has accumulated 
enough mass, the evolution of T beyond the deceleration radius 
follows the conventional decay T(R) oc R 1 appropriate for a 
parabolic flow. 

3.3. The X-Ray Flux 

The fact that the Lorentz factor remains constant for a period 
does not guarantee that the corresponding X-ray flux does 
too. Such an outcome depends also on the particular process 
responsible for the X-ray emission. For example, if the RBW 
propagates in a medium of constant density (as assumed herein), 
with the magnetic field in equipartition with the plasma (i.e., 
B — 0.4(/i()/ 1 cm _3 ) 1/,2 r G), the observed flux would increase 
with time because the source specific intensity would remain 
the same (it sweeps the same amount of electrons per unit 
time) while its solid angle (its size) increases. Therefore, the 
constant XRT emission during the plateau stage imposes certain 
restrictions on the emission process. If the observed X-ray 
emission is due to bulk Comptonization, just like the y-ray 
prompt emission, then the R~ 2 decrease of the ambient photons 
could indeed offset the oc R 2 increase of the RBW surface to 
produce a constant X-ray flux. An alternative is that the plateau 
X-ray emission is due to synchrotron emission but by a magnetic 
field that decreases with radius like B{R) oc R 1 , since the 
emissivity is proportional to B(R) 2 oc R~ 2 , thereby offsetting 
again the increase in the RBW area. Finally, if the ambient 
particle density decreases like R 2 . i.e., the RBW propagates in 




Figure 3. Lorentz factor evolution of two conical RBWs propagating in the constant density environment of Figure 1. These are distinguished from that of Figure 1 
only by their initial rest mass energies, which are Mqc 2 = 3.5 x 10 51 and 2.0 x 10 51 erg for the left and right figures, respectively. The dashed curves represent the 
evolution with radius in the adiabatic case without radiative drag. The remaining GRB parameters are given in the text. 

(A color version of this figure is available in the online journal.) 
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Figure 4. Bolometric luminosity of a RBW propagating in a constant density 
environment for a parabolic GRB outflow. The rising part is the prompt BAT 
emission, while the steep and plateau parts are the XRT afterglow, one assuming 
that the density of the photons scattered in the “mirror” decreases as R 2 . The 
time to represents the time of the onset of the burst. 

(A color version of this figure is available in the online journal.) 

a stellar wind environment, that would also lead to a constant 
flux (Shen & Matzner 2012); however, the present calculations 
for the evolution of F would have to be revised to reflect the 
different density dependence on R. 

The situation is not too different in the parabolic expansion 
case, provided that the observer “sees” only a fraction of the 
expanding RBW front. In Figure 4, we show the bolometric 
light curve of a GRB with the kinematic characteristics of 
the parabolic outflow shown in Figure 2, assuming that the 
emission past the radiative reaction slow down is due to bulk 
Comptonization, i.e., that the number of photons decreases 
with radius as R~ 2 . One has to bear in mind, however, that 
this is only a toy model light curve, based on a number 
of simplifying assumptions. However, this, along with the 
dynamics of radiation reaction feedback, set the stage for a 
more detailed future exploration of these issues. 

4. SUMMARY AND DISCUSSION 

In the present work, we have studied the evolution of the 
Loren tz factor T of a GRB RBW from its origin (To ~ I at R = 
Ro) through its acceleration (T oc R for conical outflows and 
r oc R 1 ' for parabolic outflows) to saturation (T = constant) 
and decay phases, within the context of the GRB SPM. We 
have argued within the framework of well-understood, defined 
physical processes and without the introduction of a posteriori 
assumptions that this evolution provides, among others: (1) A 
definition of what constitutes the prompt GRB phase, its broader 
spectral features, and a criterion for its termination and the onset 
of the afterglow. (2) An account of the steep decline and plateau 
or the steep decline and power law phase of the GRB afterglow, 
observed in the largest number of afterglows (Evans et al. 2009) 
and their relation to the prompt emission properties. 

Of these, (1) has been discussed in several of our previous 
publications (Kazanas et al. 2002; Mastichiadis & Kazanas 
2006, 2009), where it was shown that the SPM can provide both 


an efficient, rapid conversion of the energy stored in relativistic 
protons in the RBW of a GRB into electrons while at the same 
time producing a spectrum with // [lk ~ m e c 2 , in agreement 
with observations. In the present work, this process has been 
placed within the broader context of an evolving RBW and it is 
shown that the prompt GRB phase lasts as long as this efficient 
transfer of energy from protons into electrons is allowed by 
the kinematics of the pair production process, which depend 
crucially on T. With the reduction of F below the critical value 
set by the py —> p e + e~ reaction threshold, this transfer stops 
and the GRB luminosity decreases precipitously along with the 
value of Lpk . 

Point (2) is the novel aspect of the present work that shows that 
incorporating the radiation reaction associated with the flowing 
of the RBW through its upstream scattered radiation (or, for 
that matter, any sufficiently intense ambient radiation) in the 
evolution of F can produce a small, sharp, (over a distance 
A R < R) but important decrease of the RBW Lorentz factor 
value r. Although small, this decrease is important because it 
pushes the value of T below that of the kinematic threshold of 
the SPM (at least within the confines of the simplified treatment 
of radiation emission used herein). As a result, the transfer of 
energy from protons into electrons ceases and leads to a steep 
reduction in the GRB luminosity by roughly a factor ~m p /m e , 
consistent with the decrease in luminosity between the prompt 
and afterglow GRB luminosities. In Figure 5 (left) and (right), 
we present a sample of two such XRT light curves, namely those 
of GRB 110420A and GRB 1202 13 A, taken from Nat Butler’s 
compilation (http://butler.lab.asu.edu/swift/older.html). In these 
figures we note with the thick yellow arrow a range of ~2000 in 
flux to provide a visual estimate of the change in flux between the 
prompt and afterglow stages, which indeed is consistent with 
the estimate given here. We have also estimated by inspection the 
same ratio in a number of other bursts in the same list; we found 
some of them to be smaller and some larger. Smaller values come 
from in cases where not all protons are “burnt,” as happens to 
be the case in the numerical simulations of this process. Larger 
values of the ratio will result from a significant decrease in 
r, considering that the observed luminosity is proportional to 
~ r 4 . Such is the case in Figure 5 (left) and (right), which 
produce values T 4 ~ 10 and 100, respectively. A more detailed 
statistical analysis of this issue is currently under consideration. 

As long as the RBW radius at the point of this transition is 
smaller than its deceleration radius Ro, following this decrease 
in r (and L) due to the radiation reaction process, the internal 
energy to rest mass ratio of the RBW is sufficiently high 
to ensure that its Lorentz factor F remains constant until its 
radius reaches the value R D . At this point, F begins its more 
conventional decline. We have also argued that during this period 
of constant T. the value of the X-ray flux can remain roughly 
constant, in agreement with observations, provided that certain 
conditions are fulfilled (described in Section 3.3). Interestingly, 
synchrotron emission from a RBW propagating into a medium 
of constant density and magnetic field would yield a flux 
increasing with time and it is excluded in most cases. The effect 
of radiation reaction feedback was considered in Mastichiadis 
& Kazanas (2008) and in more detail in Mastichiadis & Kazanas 
(2009). However, the conditions of the Mastichiadis & Kazanas 
(2009) study were such that the radiation reaction transition 
radius was close to R D , just like that of Figure 3 (right), so the 
constant T section was not discernible (nonetheless, the steep 
and then more shallow decline of the 10 keV flux was apparent 
in their Figure 2 (right)). 


6 


The Astrophysical Journal, 779:16 (8pp), 2013 December 10 


Sultana, Kazanas, & Mastichiadis 



Figure 5. Flux evolution of the prompt ( Swift BAT, red points) and the afterglow ( Swift XRT, blue points) of two GRBs, namely GRB 1 10420A (left) and GRB 12021 3 A 
(right). Both GRBs exhibit the typical steep decline and plateau shape of afterglows discovered by the XRT aboard Swift. The yellow arrows span a flux range of 
approximately ~m p /m e , the value implied by the physics of the SPM. 

(A color version of this figure is available in the online journal.) 


Viewed more broadly, Figures 1 and 3 (left) and (right) 
present a set of calculations of the evolution of the Lorentz 
factor of three different RBWs with the same initial conditions 
of radius ( R t] = 10 11 cm), internal energy ( Eq = 1(P 4 erg), 
and background (uniform) density ( n = 100 cm -3 ), but with 
different values of their initial internal rest mass energy Mqc 2 , 
namely Mqc 2 = 5 x 10 51 , 3.5 x 10 51 , and 2 x 10 51 erg, 
respectively. The following values of asymptotic Lorentz factor. 
Too corresponds: 200, 286, and 500, respectively. These values 
are indeed achieved and when both thresholds of the SPM are 
satisfied, the proton energy is released to produce the main 
GRB emission and, under the force of the radiation reaction, 
the Lorentz factor is reduced sharply to T ~ 160 in all three 
cases. However, the duration of the plateau in the evolution of T 
becomes shorter with decreasing Mqc 2 . This is an important fact 
in view of the systematics between the plateau X-ray luminosity 
and its duration 7b r k (Dainotti et al. 2008, 2010). The goal of 
the present paper is not to resolve this issue but to point out 
that the considerations discussed herein, which associate R D 
with 7bk rather than with the onset of the afterglow, provide 
a novel framework within which such issues can be discussed 
and perhaps resolved. Furthermore, since emission both prior 
to and after the radiation reaction reduction in T takes place 
while r was approximately constant, it may not be surprising to 
find that the GRB properties in these two stages are correlated, 
even though one stage belongs to the prompt epoch and the 
other belongs to the afterglow GRB stage. This is indeed the 
correlation found by Sultana et al. (2012). 

One of the most challenging issues of the SPM, raised a 
number of times at conferences and indeed by the referee 
of this paper, is that of the fast variability. At <<C t, of 
the GRB prompt emission, a prominent characteristic of the 
majority of GRBs. This poses a problem, considering that 
the propagation of a blast wave in a uniform medium cannot 
produce variability times shorter than R/cT 1 . Considering 
small, spherical inhomogeneities in the swept-up medium was 
shown to be very inefficient (Sari & Piran 1997). The broadly 
accepted solution, as of this writing, is that of “internal shocks” 
(Rees & Meszaros 1994), i.e., variations in the activity of 
the GRB central engine that result in different sections of the 


outflow catching up with each other and colliding to produce the 
observed variability. However, Narayan & Kumar (2009) argued 
that this process is only moderately efficient (~1%-10%); at the 
same time, they noted that model fits to the data of specific 
GRBs produce prompt emission locations much larger than 
those inferred from the colliding shock locations (as does the 
SPM). For these reasons, these authors opted for relativistic 
turbulence to provide additional boosting to small regions within 
the broader blast wave, indicating that this notion could resolve 
the GRB fast variability issue. On the other hand, numerical 
simulations by Zrake & McFadyen (2013) have shown that 
such turbulence would quickly dissipate in the absence of 
a continuous energy input. It is possible that such an input 
is provided by the reconnection of the turbulently amplified 
magnetic field in the post-shock region, as suggested by Zhang 
& Yan (2011), who presented a detailed investigation of GRB 
variability based on the notion of magnetic field reconnection. 

One could also speculate that these processes may be possible 
to integrate within the SPM, a process that is beyond the 
scope of the present work. Clearly, a uniform medium cannot 
produce the observed variability and a collection of “blobs” in 
the ambient medium can only be efficient under the conditions 
discussed in Narayan & Kumar (2009). Nonetheless, recent PIC 
simulations (Silva et al. 2003; Nishikawa et al. 2006; Ramirez- 
Ruiz et al. 2007) have shown that relativistic shocks such as 
those from GRB RBWs are unstable to filamentation via Weibel 
(Weibel 1959) or two-stream (Buneman 1958) instabilities. If 
the column density of the RBW is 2, the splitting of the flow 
by these instabilities into, say, N filaments, will increase their 
local column density to ~/VX (there will likely be a distribution 
of columns, each becoming supercritical at different times); 
since the crucial quantity for the SPM is the local column, 
filamentation makes the conversion of RBW kinetic energy into 
radiation all the more efficient, with each such filament playing 
the role of the turbulent eddies in the model of Narayan & 
Kumar (2009). Of course, each such filament independently 
produces (via the SPM process) relativistic electrons that plow 
through the ambient radiation to produce the observed emission, 
which is perceived as a single pulse in the overall GRB prompt 
emission. Admittedly, these PIC simulations probe only scales 
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associated with their microscopic plasma physics quantities; 
however, Silva et al. (2003) do contend that these filaments 
may eventually organize into much larger scales (of the order 
of the shock width in our case). Whether these can reproduce 
the observed GRB variability is, at present, an open issue. The 
duration of the shots resulting from each such filament should 
then depend on the local electron cooling time. 

In summary, we note that the SPM, when integrated within 
the entire evolution of a GRB blast wave, provides, to the best of 
our knowledge, the first consolidation of the broader temporal 
and spectral properties of both the prompt and afterglow GRB 
stages within the framework of a single model. Considering that 
this model involves essentially no free a priori assumptions, it 
should not be expected to account for specific features of specific 
bursts; rather, it should be considered as a broader framework 
within which one attempts to account for the more specific GRB 
systematics. In this respect, we find it extremely encouraging 
that the same physics that provide for the efficient conversion 
of blast wave kinetic energy into radiation and the value(s) 
of .Epk also provide an account of the vexing XRT light curves. 
Furthermore, it makes a prediction for the bolometric luminosity 
prior to and after the radiation reaction reduction in T, indicating 
that this should be of order m p /m e . Indeed, a cursory search 
through the combined BAT-XRT light curves suggests that this 
is indeed the case; we are currently involved in providing a more 
complete statistic of the above statement, which will appear in 
a future publication. 
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